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We calculate the damping of the Bogoliubov- Anderson mode in a one-dimensional two-component 
attractive Fermi gas for arbitrary coupling strength within a quantum hydrodynamic approach. 
Using the Bethe-Ansatz solution of the ID BCS-BEC crossover problem, we derive analytic results 
for the viscosity covering the full range from a Luther-Emery liquid of weakly bound pairs to a 
Lieb-Liniger gas of strongly bound bosonic dimers. At the unitarity point, the system is a Tonks- 
Girardeau gas with a universal constant uq = 0.38 in the viscosity £ = achn for T = 0. For the 
trapped case, we calculate the Q-factor of the breathing mode and show that the damping provides 
a sensitive measure of temperature in ID Fermi gases. 



I. INTRODUCTION 

In the past few years, ultracold atoms have entered a 
new regime, where strong correlation effects appear even 
in extremely dilute gases. Prominent examples for this 
new area in atomic physics are the study of the crossover 
from a BCS-type superfluid of extended Cooper pairs to 
a BEC of strongly bound molecules [l], 0, Q or the real- 
ization of a Tonks-Girardeau gas of hard-core Bosons in 
one-dimensional atomic wires jj, H| • In the first case, the 
strong interaction regime is reached in a direct manner 
because the scattering length a near a Feshbach reso- 
nance becomes of the same order or even larger than the 
average interparticle spacing kp 1 . In the second case, 
it is the squeezing of the kinetic energy in an optical 
lattice which enhances the role of interactions [6(. A 
unique role in the context of strongly interacting ultra- 
cold gases is played by the so-called unitary Fermi gas, 
where the dimensionless interaction strength parameter 
fcpa is infinite. Thisproblem was originally discussed in 
nuclear physics @, [!]• In its simplest form, it consists 
of a two-component Fermi gas with a zero range attrac- 
tive interaction which is just about to bind a state at 
the two-particle level. Such a situation is realizable with 
cold gases at a Feshbach resonance, where the scatter- 
ing length diverges Q. Precisely at this point and for 
broad Feshbach resonances, where the range of the effec- 
tive interaction is much smaller than the mean interpar- 
ticle spacing, the full many-body problem has the bare 
Fermi-energy ep as the only energy scale. As a result, 
the complete thermodynamics is a universal function of 
the ratio ksT/cF While a quantitatively reliable 

description of the many-body problem near a Feshbach 
resonance at finite temperature is still an open problem 
[ill Il2| , the situation near zero temperature may be un- 
derstood in a straightforward manner. Indeed, at low 
temperatures, a two-component Fermi gas will be in a su- 
perfluid state, independent of the strength of the attrac- 
tive interaction. On quite general grounds therefore, the 
low lying excitations above the ground state are sound 



modes of the Bogoliubov- Anderson type, which are the 
Goldstone modes of the broken gauge symmetry in a neu- 
tral superfluid. In this regime, an effective low energy 
description is possible in terms of a quantum hydrody- 
namic approach [l3| . For a Fermi gas with a short range 
attractive interaction the associated effective field theory 
was recently discussed by Son and Wingate [lj] . Starting 
from a Lagrangian formulation of the many-body prob- 
lem, they realized that in the particular case of a unitary 
Fermi gas, there is an additional conformal symmetry 
with respect to arbitrary reparametrizations of the time. 
Remarkably, the effective field theory can be extended 
to nonequilibrium problems within the framework of lin- 
ear, irreversible thermodynamics. In particular, confor- 
mal invariance of the unitary Fermi gas applied to the 
dissipative part of the stress tensor requires that two of 
the bulk viscosity coefficients vanish As a result, no 
entropy is generated in a uniform expansion. The exten- 
sion of the effective field theory to irreversible processes 
makes evident that not only the thermodynamics but also 
dynamical properties like the kinetic coefficients are uni- 
versal at the unitarity point, a fact, first emphasized by 
Gelman, Shuryak and Zahed [16]. An example of partic- 
ular interest is the shear viscosity n which determines the 
damping of sound and collective oscillations in trapped 
gases IHj]- At unitarity, its dependence on density 
n and temperature T is fixed by dimensional arguments 
to be r\ = hna(T/fi), where fi is the chemical potential 
and a(x) a dimensionless universal function [151 ]. At zero 
temperature, in particular, r/(T = 0) = a^hn is linear in 
the density with a universal coefficient a v . Using a sim- 
ple fluctuation-dissipation type argument in the normal 
phase, a lower bound of the form a v > l/6ir has been 
derived by Gelman, Shuryak and Zahed [lg|, in analogy 
to rigorous bounds for the ratio n/hs > l/An between the 
viscosity r\ and the entropy density s in supersymmetric 
pure gauge Yang-Mills theories [H, H(| • Based on these 
results, it has been speculated that ultracold atoms near 
a Feshbach resonance are a nearly perfect liquid [H| . 

In the present work, we calculate the viscosity £ of 
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a strongly interacting Fermi gas in the whole regime of 
coupling strengths for the particular case of one dimen- 
sion, where an exact solution of the BCS-BEC crossover 
problem has recently been given using the Bethe-Ansatz 
|2ll [22. 23]. It is shown that, at T = 0, the viscosity has 
the form £ = a^hn with a coupling dependent constant 
which takes the universal value = 0.38 at the uni- 
tarity point. At finite temperature, the sound damping 
does not have a hydrodynamic form and increases like 
y/T. We determine the resulting damping of the breath- 
ing mode in a trapped gas and show that its Q-factor 
provides a sensitive measure of temperature in strongly 
interacting ID gases. 

II. SOUND DAMPING AND VISCOSITY OF A 
ID SUPERFLUID 

A. BCS-BEC Crossover in ID 

Our calculations are based on an exactly solvable 
model of the BCS-BEC crossover in one dimension pro- 
posed by Fuchs et al. [H SH and by Tokatly [2l|. 
The underlying microscopic Hamiltonian is that of the 
Gaudin-Yang model [24| 

h 2 N d 2 

i— 1 1 i<j 

of a spin 1/2 Fermi gas interacting via a short range po- 
tential giS(x). Here N is the total number of Fermions 
and m their mass. At zero temperature, the model is 
characterized by a single dimensionless coupling constant 
7 = mgi/h 2 n, where n = N/L is the ID density. For at- 
tractive interactions, the Hamiltonian (TT]) describes a so- 
called Luther-Emery liquid. Its ground state at 7 — > 0~ 
is a BCS-like state with Cooper pairs, whose size is much 
larger than the average inter-particle spacing. With in- 
creasing magnitude of 7 one reaches the strong coupling 
regime of tightly bound molecules which behave like a 
hard core Bose gas as 7 — > — 00. As shown by Gi- 
rardeau, the hard core Bose gas in one dimension is 
equivalent - for densities diagonal in real space - to a 
gas of non-interacting Fermions [25j . Within a strictly 
one-dimensional model, the BEC-limit of strongly bound 
pairs is thus a Tonks- Girardeau gas. Now in practice, the 
atoms are trapped in a harmonic waveguide with radial 
frequency uj±/2tt. The associated transverse oscillator 
length a± = \J h/muj± then defines an additional length, 
not present in the Gaudin-Yang model |T]). As shown by 
Bergeman et al. [26| , the exact solution of the scattering 
problem for two particles in such a waveguide, interacting 
with a 3D pseudo-potential with scattering length a, al- 
ways exhibits a two-body bound state, whatever the sign 
and magnitude of the scattering length a. It appears at 
an energy huj±_ — e&, which is below the continuum thresh- 
old at Huj± of the transverse confining potential. Apart 
from this bound state, all the scattering properties can 



be described by an effective ID delta potential g^ a S(x) 
for atom-atom interactions with strength [27j 

gl a = 2hw ± a{\ - Aa/aA_)~ X . (2) 

As naively expected, an attractive 3D interaction a < 
implies a negative value of g^ a . The associated binding 
energy tb = m(g1 a ) 2 /Ah 2 in the ID delta potential co- 
incides with the exact value in the weak confinement 
limit \a\ -C a±. Remarkably, the strength gf a of the ID 
pseudo-potential remains finite at a Feshbach resonance 
where a = ±00. The corresponding exact value of the 
binding energy is lb ~ 0.6ftw± [2l|,|26j]. Entering the pos- 
itive side a > of the Feshbach resonance, the vanishing 
of the denominator in ^ at aj_/a = A ~ 1.0326 [281 ] 
leads to a confinement induced resonance (CIR), where 
g1°" jumps from — 00 to +00. The exact bound state at 
this point has binding energy = 2huj± and a spatial 
extension along the x-axis, which is of the order of the 
transverse oscillator length a±. With decreasing values 
a < a± of the 3D scattering length, tb increases mono- 
tonically beyond 2huj± and finally approaches the stan- 
dard 3D result h, —> h 2 /ma 2 in the weak confinement 
limit a <C a± [U, [2||. Since fkj± > £f in the limit 
of a singly occupied transverse channel, the true bound 
state energy q, is the largest energy scale in the prob- 
lem in the regime after the CIR where gf a > 0. In this 
regime, the appropriate degrees of freedom are no longer 
the single atoms but instead are strongly bound Fermion 
pairs, which are essentially unbreakable. An exact so- 
lution of the four-body problem in a quasi ID geome- 
try with tight harmonic confinement shows, that these 
dimers have a repulsive interaction in the regime beyond 
the CIR HI. The related constant gf d > in the effec- 
tive dimer-dimer interaction gf d 5(x) can be calculated as 
a function of the 3D scattering length 29] . It approaches 
gf d — > 2huj± ■ 0.6 a — > in the weak confinement limit, 
where the dimer-dimer scattering length a dd ss 0.6 a is 
identical with the one in free space [301 ]. Sufficiently far 
from the CIR, one thus recovers a weakly interacting gas 
of dimers. 

At the many-body level, the situation after the CIR 
is described by a Lieb-Liniger model (3~H of repul- 
sive Bosons. Its dimensionless coupling constant 7 = 
mg dd /h 2 n is now positive and vanishes in the weak con- 
finement limit. It diverges at a value of the 3D scat- 
tering length a of order a±. Now although the diver- 
gence of gf d does not exactly coincide with that of gf a 
at the CIR [29( , the range of inverse dimensionless cou- 
pling constants where this mismatch appears is of or- 
der 1/7 « na± [23j]. It is thus negligible in the rele- 
vant low density limit (naj_) 2 <C 1. Indeed, at a fixed 
density n, the quasi ID condition Huj± tp that only 
the lowest transverse mode is occupied is equivalent to 
(naj_) 2 « 1. In the limit na± — * 0, there is thus a 
continuous evolution from the Gaudin-Yang model of at- 
tractive Fermions to the Lieb-Liniger model of repul- 
sive Bosons which completely describes the BCS-BEC 
crossover in one dimension [U [22j ■ The associated spec- 
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trum of elementary excitations is straightforward to un- 
derstand: in the BCS limit I/7 — * —00, the system con- 
sists of weakly bound Cooper pairs. Their breaking is 
associated with a finite excitation gap and the corre- 
sponding spectrum exhibits a relativistic dispersion re- 
lation e s (k) = \J (A/2h) 2 + (v s k) 2 similar to the stan- 
dard quasiparticle spectrum of the BCS theory. The as- 
sociated energy gap A and the spin velocity v s > vf 
increase monotonically with I/7, both diverging in the 
strong coupling limit 1/7 = at the CIR [22|]. In ad- 
dition, there are gapless density fluctuations describing 
the Bogoliubov-Anderson mode of a neutral supcrfluid. 
These excitations exist for arbitrary coupling, both be- 
fore and after the CIR. Their spectrum is e(k) — v c \k\ at 
low momenta, with a (zero) sound velocity v c [32J , which 
monotonically decreases from the ideal Fermi gas value 
v c = Vf at 7 — > 0~ to the weak coupling BEC result 
v c = Vjvp/ir as 7 -> 0+ [H. A t the CIR, 1/7 = 0, the 
system is a Tonks- Girardeau gas [I^] of tightly bound 
dimers. The value u c (l/7 = 0) = Uf/2 simply reflects 
the fact that the unitary Fermi gas in ID is a hard core 
Bose gas which - in turn - behaves like an ideal gas of 
spinless Fermions at half the original density. The uni- 
versal parameter /3, which follows from v c = Upy! + (3 
[l(J, thus has the exact value = — 3/4 in one dimension. 



B. Quantum Hydrodynamic Theory 

In order to calculate the damping of long-wavelength 
phonons, we use a one-dimensional version of quan- 
tum hydrodynamics (QHD). Restricting the attention to 
the gapless Bogoliubov-Anderson mode in the superfluid 
regime, the QHD Hamiltonian has the same form on both 
the Fermionic (before the CIR) and the Bosonic (after the 
CIR) side of the ID BCS-BEC crossover. In a harmonic 
approximation, which is valid at low energies, the sound 
mode is described by the quadratic Hamiltonian 

Ho= v i r<ixP(M 2 + -n 2 ] (3) 

2 Jo lv c Po > 

where the conjugate fields f{x) and Tl(x) describe phase 
and density fluctuations respectively. The only input pa- 
rameters are the equilibrium mass-density pa — mn and 
the sound velocity v c . From the Bethe Ansatz, the veloc- 
ity v c is known as a function of the dimensionless inverse 
coupling constant 1/7, which ranges between 1 Ay = — 00 
in the BCS- and I/7 = +00 in the BEC-limit To de- 
termine the damping due to the interaction of phonons, 
the energy functional of a one-dimensional quantum liq- 
uid needs to be expanded beyond quadratic order in 
the fields <p(x) and H(x). It is a crucial advantage of 
the QHD approach, that the coefficients of the leading 
nonlinear terms are completely determined by thermo- 
dynamic quantities fl3j . Specifically, the lowest (third) 
order terms give rise to a contribution H, m t to the total 



Hamiltonian of the form 

H in t = ^J dx {(d x <p)n(d x <p) + (d x ^fn 

™ 2+ i(f) n3 } (4) 

The quadratic Hamiltonian (J3]) is diagonalized by the 
standard mode expansions 




where b q and b q denote the usual bosonic creation and 
annihilation operators respectively and Vp = nhn/2m 
is the Fermi velocity of the noninteracting gas. After 
inserting the mode expansions in we obtain 

H int = ^2 -j= V(qi,q 2 ,q3) { b qi b q2 b q3 S(qi + q 2 + q 3 ) 

91 192,93 

+ b qi b q2 b \j{ll + 12 ~ 33) + b qi b l 2 b q 3 S (ll ~ Q2 + (ft) 

+ ML 6 l%i - 22 - ?s) + fr-c } (5) 
with the vertex 



V(qi,q2,qa) 




The momentum dependence of the vertex makes pertur- 
bation theory applicable for long wavelength phonons. 
More precisely, it applies as long as the imaginary part 
of their energy is much smaller than the real part. As will 
be shown below, this requires the phonon wavelength A to 
be much larger than the interparticle spacing on the BCS 
side of the crossover, while in the BEC limit the more 
resctrictive condition An 3> (I/7) 1 / 4 ^ 1 is required. It 
should be remarked that an approach to evaluate the non- 
linear terms in the Hamiltonian via Bosonization leads to 
additional coupling terms between phonons and spin ex- 
citations on the BCS side of the crossover which are not 
accounted for in the QHD approach. As will be shown 
in the Appendix, however, these terms do not contribute 
to the phonon damping at long wavelengths. 

We evaluate the damping by calculating the imaginary 
part of the phonon self-energy S(fc, uj) which is defined as 
the analytic continuation of the corresponding self-energy 
Yi th in the exact thermodynamic Green function 

Q k,iw n i Un - Vc \k\ -T, th (k,iu) n ) ' 
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Here ui n = 2irn/[3 with n S Z are the standard Bosonic 
Matsubara frequencies. (We set h = kg = 1 from now 
on, except in final results) 

The main contribution to the damping rate comes from 
the three self-energy diagrams shown below, correspond- 
ing to spontaneous decay, absorption of a phonon and 
three wave annihilation. 

k — q 




q — k 




Here G R and G A are the usual retarded and advanced 
Green functions respectively. The combinatorial factor 
18 arises from the three possible ways of choosing the 
creation/annihilation operators for the initial and final 
phonon and two possibilities of pairing the phonons in 
between. As was already pointed out by Andreev and 
Samokhin, the fact that for linearly dispersing phonons in 
ID, momentum- and energy conservation are simultane- 
ously satisfied, requires to go beyond second order pertur- 
bation theory which would give an infinite damping rate. 
Following the approach of Andreev [33| , we calculate the 
self-energy by using the fact that T, R (k,uj) <C e k := u c |fc| 
at long wavelengths. The precise condition on k for which 
this holds, has to be determined afterwards and will be 
discussed below. Since we are interested in the quasipar- 
ticle pole of G R (k,u) = — E fl (fe, ui)) we can use 
the approximation 

ui = e k + E fl (fc, w) w £ fe + Y, R {k, e k ) =: e k + E* 



-k — q 




Diagrams of this type have been considered before by An- 
dreev [2H who studied the sound absorption in ID Bose 
liquids for T > and by Samokhin 34] in the context 
of the damping of zero sound in a ID liquid of repulsive 
Fermions. 

After taking the limit L — > 00 and analytic continu- 
ation, the retarded self energy H R (k,oj) is given by the 
sum of the three diagrams 

E fl (fc, u) = Ef (k, oj) + 2 E£(fc, oj) + Ef (fc, «) (7) 

with 

Ef(fc,c) = -18 J coth(^) V 2 (k,q,k-q) 

x { - G R (q, ft + cj) Im G R (k - q, -ft) 
+ G R (k-q,LU- fl) Im G R (q,fl)} 

X R (k,Lo) = -18 J g^g COth(^) V*(k,q,q-k) 
x ^G R (q, fl + u>) Im G R (q - k, ft) 
+ G A (q -k,fl-w)Im G R (q, ft)\ 



Ef(fc,c) = -18|ggcoth(^) V\k,q, 
x I - G A {q, ft - ui) Im G R (-k - q, 
+ G A (~k -q,-u- ft) Im G R {q, Q)| . 



-k-q) 
-ft) 



leading to 



G R (k,uj) 



-e k -^ R 



(8) 



The damping rate of phonons with wavevector k > 
can now be determined from the imaginary part T k := 
Im 'E (k, £fc) = Im Yj R of the (on-shell) self energy. 

Starting with the case of zero temperature T = 0, the 
only contribution to E^ comes from spontaneus decay 
(first diagram). Applying the approximation |[8]) in the 
integrand and doing the ft- integration, one ends up with 



dq 
2^ 



+ G R (k-q,e k 



V 2 (k : q,k-q)[G R (q, 

>}■ 



£fc — £fc- 



yR 



The major contribution to the integral comes from < 
q < k where e k — e k - q 
equation 



e q = 0. Thus we arrive at the 



Ei = - 



-ITVr 



32 m 2 i>F 
k 



3+ 4 _d_ 

v 2 dvp \vp 



f p-kq(k-q) 

Jo 27r 



yR 



for the retarded self energy. It is solved with a purely 
imaginary Ansatz E^ = — i/it? 2 where 



Am 



7T V c 

2v F 



1 



/i(o = 2W3 + 



dx- 



v 2 dvp \v F 



x(l — x) 



2tt Jo ~~ (1 - x) a + x a 
At zero temperature, the resulting damping rate 

v 2 d 



n 



8m y vf V 4 



v 2 dvi 



k 2 (9) 
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of the Bogoliubov- Anderson mode in one dimension is 
therefore quadratic in the wavevector. Formally, this is 
precisely the behaviour of a hydrodynamic mode. It al- 
lows to define a zero temperature viscosity £ by the rela- 
tion 



form 



uok = v c k - 



C 



2mn 



k 



(10) 



which is completely analogous to sound damping in three 
dimensions. In that case, £ is replaced by the combina- 
tion £2 + 477/3 involving one of the superfluid bulk vis- 
cosities C2 and the shear viscosity r] [IB, [35| and one has 
£2 = at unitarity [15]. From the result ((9]), we see that 
the viscosity at zero temperature has the form 



with a constant 




Vc_ 



d 
dvp 



(11) 



(12) 



A plot of is given in figure [TJ where the exact result 
for v c /vf from the Bethe-Ansatz solution was used. Evi- 
dently, the dimensionless viscosity coefficient depends 
on the inverse coupling constant 1/7 of the BCS-BEC 
crossover. It is thus in general dependent on the par- 
ticle density n. At the unitarity point, however, where 
1/7 = 0, this dependence vanishes and takes the uni- 
versal value 



0.38 



which is just l/v2 of the value 0^(7 = 0~) = 0.54 at- 
tained in the weak coupling limit of the ID noninteracting 
Fermi gas. Concerning the range of applicability of the 
perturbative calculation, it is obvious that the approxi- 
mation leading to © is satisfied for small wavenumbers 
k -C mnVcl Based on the explicit result for £ this con- 
dition translates into phonon wavelengths much larger 
than the mean interparticle spacing on the BCS side of 
the crossover, including the unitarity point. On the BEC 
side, the ratio u c /C vanishes like y/v^ ~ 7 1 / 4 . The phonon 
wavelenghts have thus to obey the more restrictive con- 
dition An ^> (I/7) 1 / 4 ^s> 1 mentioned above. It is inter- 
esting to note, that this condition is less restrictive than 
the requirement A ^> £1 for the validity of the linear dis- 
persion relation ujk = v c \k\, where £1 = n _1 (l/7) 1 / 2 is 
the ID healing length. 

We now turn to the situation at finite temperature T > 
0, where the long wavelength phonons, for which T 3> 
behave classically. The thermal factor coth(/3Q/2) may 
therefore be replaced by its classical limit 2 /(/3ft). At 
T / 0, the second diagram representing the absorption 
of another phonon also contributes to the damping. An 
explicit calculation along the lines performed at T = 
gives a phonon damping rate at finite temperature of the 



1 k 



" fcsT / 2 (3/2)h ' 



Am V 2 hvp 



I d 



v?. dvi 



v F 



fc 3 / 2 . 
(13) 

Here /2(a = 3/2) w 0.6221 is a numerical coefficient de- 
fined by the integral 



h(a) := 



1 dx 



1 



2w (1 - x) a + x a J x 2tt (x - l) a + x a 



dx 



1 



As already noted by Andreev [33)], the damping ~ fc 3 / 2 
for T > is not of the standard hydrodynamic form, in 
contrast to the behaviour at zero temperature. The quite 
different results can be understood from the fact that at 
any finite temperature, the quasi long-range superfluid 
order present at T = is destroyed by phase fluctuations 
on a characteristic length scale £t := hv c /ksT. Depend- 
ing on the ratio y — k^T between this length scale and 
the phonon wavelength, the behaviour is either essen- 
tially superfluid for y > 1 or normal for y <C 1. Similar 
to the formulation used in dynamical scaling laws near 
critical points [H, [37j , the crossover between the two dif- 
ferent types of behaviour may be described by an Ansatz 
of the form 



hk 2 

T k = ^^ T k). 
2m 



(14) 



The associated crossover function has the limiting be- 
haviour 



3.70 q c 



for y —> 00 (T — > 0) 



for y < 1 (T > e k ) 



with the parameter defined in Eq. (fl2|) . It should 
be pointed out, that the dependence of the damping rate 
on temperature is a simple power law ~ T 1 / 2 only to the 
extent that the temperature dependence of the velocity v c 
itself can be neglected. Moreover, note that for nonzero 
temperature, the damping (|13[) remains finite in the BEC 
limit I/7 — > 00 in contrast to the T — case. 



C. Harmonically trapped gas 

Finally we extend our results on damping in a homoge- 
neous gas to the experimentally accessible case of a har- 
monically trapped system, using essentially a local den- 
sity approximation. Our main interest is to calculate the 
damping of the so-called breathing modes, which have al- 
ready been measured in ID Bose gases in a regime near 
the Tonks-Girardeau limit 3S|- Assuming a standard 
type of viscous damping in a classical fluid, the damping 
rate in a ID inhomogeneous case is given by [391 ] 



{-&mech)t 




Jdz az){{d z v?) t 


2(E) t 




2m J dz n(z)(v 2 )t 



(15) 
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FIG. 1: Viscosity-parameter as a function of the inverse 
coupling constant j~ . 



where (.)t denotes the time average, z is the spatial co- 
ordinate and the last equation holds for harmonically os- 
cillating perturbations, where (E) t = 2(Eki n )t- 

Breathing modes in a harmonic trap are characterized 
by a velocity profile of the form v(z, t) = const, z e~ tuJBt . 
Since the damping of the Bogoliubov- Anderson mode at 
zero temperature has precisely the form of a standard 
viscous fluid one obtains using Eq. (fTTj) 



1 B 



(16) 



2m (z 2 ) 

where the brackets denote the spatial average defined by 



(/(*)> 



1 

N 



dz n(z)f(z) 



and N is the total number of particles. Since the con- 
stant depends on density except at the unitarity point, 
the damping also involves a spatial average (a^). A plot 
of the coupling dependent Q-factor Q = °f the 

breathing mode at zero temperature is given in figure [5] 
together with the ratio of its frequency in units of the 
axial trap frequency lo z . The required density profiles 
were calculated numerically using a local density approx- 
imation and the exact results for the chemical potential 
from the Bcthc-Ansatz solution. As shown by Menotti 
and Stringari [4(| , the density profiles also determine the 
frequency from 



din (z 2 



diu 2 



The appearance of a maximum in the breathing fre- 
quency just before the confinement induced resonance 
may qualitatively be understood from the fact that 
around this point the nature of the pairing changes from 
overlapping, correlated pairs to individual molecules. In- 
deed, the size of a molecule is of the order of the average 
inter-particle distance n _1 for inverse coupling constants 
I/7 ~ —0.5 Remarkably, this is also close to the 



FIG. 2: Frequency lob and quality factor Q — lob /To of 
breathing modes at T — as a function of the 3D scatter- 
ing length a. The dashed line indicates the CIR (unitarity 
point). The plot is for u± — 5Nui z , where Wj_ and uj z are the 
radial and axial trapping frequency, N denotes the number of 
particles in the trap. 



point, where the effective dimer-dimer interaction gf d di- 
verges [29] and where the dimensionless viscosity coef- 
ficient (|12p exhibits a small maximum as shown in Fig. 

m 

For T > the situation is more complicated, because 
the k -dependence of the damping rate in (|13[) implies 
a non-hydrodynamic behaviour. In order to account for 
the specific k -dependence in a classical calculation, we 
modify the stress-'tensor' by introducing an effectively 
velocity dependent viscosity. The form of £ to be used 
in (fTS"]) which leads to the result (JT4J) for the damping in 
the homogeneous system is 



( = hn $ £t 



<M> 



(17) 



with the function $ as defined in JR 

Using (jTTJ) for the inhomogeneous system, the result 
for the damping of breathing modes can be expressed as 



h 

2m 



r B = — ($ 



6 



1 

(?) 



737 (18) 



In the particular case of the Tonks- Girardeau limit de- 
scribing the unitary ID Fermi gas at the CIR the damp- 
ing of breathing modes at T = and at finite tempera- 
tures T 3> huj z is given by 



F° 

1 B 



1 B 



1.5 



4.1 



N 

Uz_ ... 

N\'hw 



k R T 



(19) 
(20) 



A plot is given in figure [3J Note that in the Tonks- 
Girardeau limit, the constraint £r <C (\z\) simply trans- 
lates into ksT 3> hu z . The zero temperature result 



7 



r T /((D z /N) 



5 10 15 20 25 30 35 40 45 50 



FIG. 3: Damping of breathing modes at the unitarity point 
as a function of T. The damping at T = is shown for 
comparison. 



for the damping is thus only valid in the experimentally 
hardly accessible regime T <C hu z , while for realistic tem- 
peratures the damping is expected to increase like VT. 



haves like ~ k 3 / 2 and thus is not of a hydrodynamic 
form. The resulting damping of the breathing mode in a 
trapped gas has been calculated within a simple model, 
which accounts for the inhomogeneity in the case of a 
nonstandard damping. In particular, it has been shown 
that in the experimentally relevant regime T 3> hcu z , the 
damping increases like y/T, thus providing a sensitive 
measure of temperature in strongly interacting ID gases. 
Experimentally, an attractive Fermi gas near a Feshbach- 
and confinement induced resonance has been realized by 
Moritz et al. [4l|. Since the typical temperatures in 
this gas were of order T » 0.2 Tp with Tp N ■ fuv z 
and typical particle numbers are N w 100, the condi- 
tion T 3> fouJz is realized. It would be quite interesting 
therefore, to study the temperature dependence of the 
breathing mode Q-factor similar to the measurements 
performed in ID Bose gases (38|. In this context, it is 
interesting to note that for the Tonks-Girardeau gas, ex- 
act results for the dynamics have been derived at zero 
temperature by Minguzzi and Gangardt [43]. In partic- 
ular they imply zero damping of the breathing mode at 
the unitarity point, i.e. an infinite Q-factor. From our 
present results, the Q-factor is infinite only in the limit 
N — > oo but not for the finite and typically small values 
N w 50 — 100 realized experimentally. This point needs 
to be studied further. 



III. SUMMARY 

In summary, we have calculated the damping of the 
Bogoliubov- Anderson mode for an attractively interact- 
ing ID Fermi gas in the whole regime between the BCS 
and the BEC limit. At zero temperature, the damping 
is of a hydrodynamic form with a viscosity £ = a^hn. 
The associated constant a( is a smooth function along 
the crossover from a BCS type superfluid to a BEC of 
strongly bound pairs of Fcrmions, with a universal value 
= 0.38 at the unitarity point. It is remarkable, that 
a rough analysis [l6| of the experiments by Bartenstein 
et al. [HI gives a value of 0.3 for the universal viscosity 
coefficient of the 3D unitary Fermi gas at the lowest at- 
tainable temperatures. However, it is obvious that a com- 
parison between this and our result in ID is not meaning- 
ful. Nevertheless, the fact that the Bogoliubov-Anderson 
mode spectrum and velocity are hardly different between 
the one and the three dimensional case, suggests that the 
viscosity in 3D exhibits a similar dependence on the in- 
verse coupling constant 1/7= \/{kpa). The unitarity 
point would then define a minimal value of the viscosity 
on the Fermi side of the crossover, yet lower viscosities 
are expected by going further into the BEC regime. It is 
a peculiar property of the ID BCS-BEC crossover prob- 
lem, that the boundary between Fermionic and Bosonic 
behaviour is sharp and defined by the confinement in- 
duced resonance. A similar sharp separation, however, 
does not exist in three dimensions. 

At finite temperature, the sound damping in ID be- 



APPENDIX A: QHD VS. BOSONIZATION 

As mentioned earlier, the QHD approach does not give 
rise to interaction terms between spin and charge excita- 
tions. We now use Bosonization techniques to construct 
Hint on the BCS side of the crossover, where spin ex- 
citations are present, and show that these terms do not 
contribute to the damping rate. 

In order to obtain damping in the Tomonaga-Luttinger 
model, one needs to to incorporate the quadratic disper- 
sion relation in the Hamiltonian, leading to third order 
terms in the fields as was already shown by Haldane for 
a spinless gas of Fermions [44| . Since we are dealing with 
spin 1/2 particles, third order terms arise which couple 
charge and spin excitations. Bosonizing the kinetic en- 
ergy term via point splitting we find 

Hint = 7T~ ^] ( dx^d x i)ld x il> a + d^dxtpA 

= ~^/f l tefavM + n c(^c) 

+ n c (a^ c )n c + {d xVc f + 3n c [(a^ s )n s 
+ n s (d x <p s )] + 3(d x <p c ) [n 2 s + (d x <p s ) 2 ] } 

where we used the same notation as 45]. A similar re- 
sult was derived previously by Fedichev et al. in the 
context of spin-charge separation in ID repulsive Fermi 
gases 46]. To avoid confusion, it should be mentioned 
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that here d x (p plays the role of density fluctuations and 
II describes phase fluctuations in contrast to The 
charge-charge-interaction Hamiltonian given above is es- 
sentially the same as the one in Q with one subtle differ- 
ence: in the QHD approach the fl 3 term has a prefactor 
eT~^( v c/Po) which is absent in the bosonized counterpart. 
This factor is important because it prevents and thus 
the damping from going to infinity in the BEC limit. 
Since the crossover from the BCS to the BEC regime is 
continuous, the damping must also change continuously 
when one crosses the point of unitarity. This argument 
leads us to include this prefactor also in Hi nt on the BCS 
side and thus use dU) in the whole crossover regime. 

The contribution to the phonon damping rate arising 
from interaction with spin excitations in second order 
perturbation theory can be calculated from the following 
diagrams: 



k — q 
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